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Abstract 

By studying the nonperturbative contribution to the Dokshitzer-Gribov-Lipatov- 
Altarelh-Parisi and Gribov-Levin-Ryskin equation, it is found that (i) the nonper- 
turbative contribution suppresses the evolution rate at the low Q^, small-x region; 
(ii) the nonperturbative contribution weakens the shadowing effect. The method 
in this paper suggests a smooth transition from the low ("soft" ), where non- 
perturbative contribution dominates, to the large ("hard" ) region, where the 
perturbative contribution dominates and the nonperturbative contribution can be 
neglected. 

PACS numbers: 12.38 Aw, 13.60.Hb 



*Permanent Address: Department of Applied Physics, Shanghai Jiaotong University, Minhang Cam- 
pus, Shanghai 200240, P.R. China. E-mail: hebobo@public6.sta.net.cn 



The properties of parton distribution at small-x region ( x is the value of the Bjorken 
variable ) have recently been the important subject [1-5]. Recent measurements of the 
structure functions for the deep-inelastic ep scattering at HERA discovered their dramatic 
rise as x decreases from 10~^ to 10"^ [6,7]. The predictions of Gliick, Reya and Vogt (GRV) 
[8] by using the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi ( DGLAP ) evolution equation 
[9] at very low ( is the negative of the square of the four-momentum transferred 
by the lepton to the nucleon ) are in broad agreement with this result. However, the 
GRV model fails in low quantitatively, that is to say, the evolution rate is faster than 
that of the experiments. Q It then becomes a challenging problem how to determine 
the structure function at the low region. Another important question is whether the 
shadowing effect, which the Gribov-Levin-Ryskin (GLR) equation [10] describes, can be 
observed by the current experiments at HERA. 

The purpose of this letter is to study the DGLAP and GLR equation in low by 
considering the nonperturbative contribution. It will be showed that (i) the nonpertur- 
bative contribution suppresses the evolution rate at the low Q^, small-x region; (ii) the 
nonperturbative contribution weakens the shadowing effect. So the nonperturbative con- 
tribution is very important in the low region. 



The DGLAP equation for the gluon distribution at small-x region in the DLLA is 
given by 



d'^xg{x,Q'^) 3as{Q 
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By considering the shadowing effect, the DGLAP equation can be modified in the 
form ( called GLR equation): 

d^xg{x,Q^) 3a,(g2) 81a|(g2) 

-xg{x,Q ) - ————[xg{x,Q )] (2) 
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^ Ref. [7] shows that the value of F2 given by GRV model is lower than that of the experiments at 
Q2 = OAGeV^ while it is higher at = 6.5GeV^. 



In the previous studies people used the perturbative QCD effective couphng ( the 
leading order coupling) agiQ"^) to study the equation, where 
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By using the formula (3), the Eq.(2) can be cast in the form: 



dydtG{y,t) = cG{y,t) - AexpH - exp(t)]G2(y, t) (4) 

where y = ln(l/x),t = ln[ln(QVA^^^)], G{y,t) = xg{x,Q^), c = 12/(11 - 2nf/3) with 
Uf the number of quark flavors and A = dTc^c^ /IGR'^Aqq^. The equation (4) has been 
studied by many authors [11,12]. Eskola, Qiu and Wang [13] applied the eq.(4) to study 
the shadowing in heavy nuclei. 

However, it should be noted that in the large the formula (3) derived from pertur- 
bative QCD is a good approximation while in the low nonperturbative contribution 
should be included. To avoid the ghost-pole problem in the behavior of a running cou- 
pling, Shirkov and Solovtsov [14] obtained the QCD running coupling in the leading order 
as: 
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The second term on the right-hand side of Eq. (5) is clearly the nonperturbative con- 
tribution. It is noted that aan{Q^ = OGeV^) = (^^^^2nf) depends only on group factors and 
does not depend on the Aqcd- It can be found that in the large the running coupling 
ctan(Q^) is dominated by perturbative contribution and the nonperturbative contribution 
can be neglected while in the low the nonperturbative contribution is very notable. 

By applying the formula (5) to DGLAP equation 0, the Eq. (1) can be written as: 

dydMy, t) = cG{y, t) [1 + "^P^^) ] (6) 

1 — exp(e*) 



^The reason that the formula (5) replaces forniula(3) in DGLAP equation will be given in the ending 
of the text. 



Adopting the semi-classical approximation [4,12], which amounts to keeping only the 
first order derivatives of the function z{y,t) — log[G{y,t)], Eq.(6) is rewritten as: 

dyz{y, t)d,z{y, t) = c[l + ^Z^^].. ] (7) 

Eq. (7) can be solved by using the method of characteristic [4,12]. Let 

p = dtz{y,t), q = dyz{y,t) (8) 

Eq. (7) can be written as the following general form: 

Fip,q,t,y,z)^0 (9) 

The characteristic differential equations of Eq. (9) have the following form: 

dr ~ ^' dr dr -P^P + ^^<i^ 

where r is the "inner" time. 

F F F^ F F K 

Fp^q, F, = p, Fy = Q, F, = 0. 

Ft = -ce*[l - exp(e*) + exp(e* + t)]/[l - exp(e*)]' (11) 

Same as Ref. [12], the initial conditions to solve Eqs. (10) are: 

to = log[log(Q^/A|c^)] = log[log(4Gel-VA|ci.)], 

yo = Iog(l/Xo) = log(l00),po = c/5bare, ^0 = log(3.38), 

^0 = + ^-TT;syPo^ ^QCD = 0.2GeV, d^^e = 0-5. (12) 

1 — exp(e*o) 



For clear comparison between the results of Eqs. (10) with those of DGLAP equation 
not including the nonperturbative contribution , we solve the Eqs. (10) numerically by 



adopting the Runge-Kutta methods [12]. Fig. 1 shows the evolution path (y,t) corre- 
sponding to the dashed line by solving the Eqs. (10) compared with the evolution path 
(y,t) corresponding to the solid line by solving the characteristic differential equations of 
DGLAP not including nonperturbative contribution. 

Like the DGLAP evolution equation, by using the formula (5), the GLR equation can 
be cast in the form: 

dyd,G{y,t)=cG{y,t)[l+ ^ ^^P^^) ,J -AexpH-exp(t)][l+ ^^^^^^ J 'G'{y,t) (13) 

1 — exp(e^j 1 — exp^e^^j 

The semi-classical approximation of the Eq. (13) 

dyz{y,mz{y,t) = c[l + "^P^^^ ] - Aexp[-t - exp(i) + z][l + "^^^^j (14) 

1 — exp^e'^j 1 — exp^e'^j 

By using the method of characteristic, the solution of Eq. (14) are showed in Fig. 2. 



To conclude, recent experiments at HERA have supplied much information about the 
nucleon structure at both large and low Q^. In large Q^, the DGLAP equation derived 
from perturbative QCD can describe the behavior of parton distribution. The challenging 
problem is how to make a unified treatment on nucleon structure at both large and low 
Q^. This letter proposes a way to meet the requirement. It is well known that the parton 
distribution includes both perturbative QCD and nonperturbative QCD effects. The 
input distribution reflects the nonperturbative QCD, and the DGLAP equation itself is 
the perturbative QCD. So it seems that the DGLAP equation describing the perturbative 
QCD effect and the input distribution describing the nonperturbative QCD effect together 
can give the comprehensive description to the parton distribution. However, it can be 
seen clearly that the input distribution does not include all nonperturbative effects, that 
is to say, some nonperturbative effects are reflected through the running coupling. By 
considering the nonperturbative effects in running coupling, it is a natural way to apply 
the DGLAP equation in the low region. So the evolution equation itself includes both 
perturbative and nonperturbative effects. It should be noted that this way is a work 
ansatz, which can not be derived from the theory. In viewing the Fig. 1, it can be found 
that the nonperturbative contribution to DGLAP is very notable. 



Although the predictions of GRV model by applying the DGLAP evolution equa- 
tion at very low are in broad agreement with HERA experiments, the evolution rate 
which results from the model is faster than that of the experiments. By considering 
the nonperturbative contribution to DGLAP, the discrepancy between GRV model and 
the experiments can be explained naturally. From analysing the DGLAP equation, it 
can be found that the running coupling determines the evolution rate, which becomes 
slow by considering the nonperturbative contribution, especially in very low such as 

Recently, one of the important questions is whether the shadowing effect can be ob- 
served by the current experiments at HERA. Some people say "can" such as Shabelski and 
Treleani [15] while other people say "cannot" such as Golec-Biernat, Krasny and Riess 
[16]. Ayala, Gay Ducati and Levin [17] argue that the shadowing effect is large in the 
gluon distribution but small in ^2(0;, Q^). Like the DGLAP equation, the GLR equation 
can be treated by the same method. In this paper, a firm conclusion about this question 
does not be made. Nevertheless, in viewing the Fig.2, it can be concluded that shadowing 
effect in the GLR equation, which is modified by the nonperturbative contribution, is not 
so notable as what has been studied in the case when gluons concentrate in "hot-spots" 
within proton {R = 2GeV~^). By analysing the value of as{Q^) and aan{Q^), the simple 
explanation of this result is that the linear term in the GLR equation is proportional 
to as{Q^) or aan{Q^) while the nonlinear term is proportional to a^Q"^) or Q!^„((5^)) so 
the shadowing effect is weakened due to q;s((5^) > aaniQ'^), especially at low Q'^ region. 
The result means that the nonperturbative contribution weakens the shadowing effect. 
Some people [13] discussed the nuclear shadowing by applying the GLR evolution equa- 
tion without considering the nonperturbative contribution. It might be interesting to 
restudy the nuclear shadowing by applying the GLR equation (14), which is modified by 
the nonperturbative contribution. 

In this paper, the DGLAP and GLR equation are solved by applying the method 
of characteristic. From the initial conditions (12), it can be found that the start point 



is = iGeV^, which is not very low, because in very low the evolution equation 
might be too complicated to treat and the semi-classical approximation is not a good 
approximation. However, the conclusions shown in this paper can be deduced easily to 
the very low region, where the nonperturbative contribution becomes more dominant 

because the difference between the formula (3) and formula (5) as(Q^) — aan(Q^) which 
is (0.5 - 0.4) = 0.1 at g2 = 0.65GeV^ is much more notable than that (0.3 - 0.29) = 0.01 
at g2 ^ 4GeV^. 

In summary, it is believed that QCD, which includes the perturbative part and the 
nonperturbative part, is a complete theory to describe all strong interaction experiments. 
Nevertheless, as the fundamental dynamical model, the DGLAP or GLR equation itself 
only reflects the perturbative effect. Unfortunately, almost all strong interaction experi- 
ments such as the deep-inelastic ep scattering at HERA involve both perturbative effect 
and nonperturbative effect. The purpose of the paper is to develop a fundamental dy- 
namical model which itself includes nonperturbative effect. By studying the model, we 
make the conclusions: (i) the nonperturbative contribution suppresses the evolution rate 
at the low Q^, small-x region; (ii) the nonperturbative contribution weakens the shad- 
owing effect. Those conclusions are helpful to explain the recent experiments in low 
region. 

If the results of this paper are compared with the recent HERA data in detail, the 
quark distribution must also be discussed in the low region, but it should be noted that 
the recent HERA data is available at a few isolated values of averaged x and Q^, especially 
in very low Q^, so how to analyse the HERA experiments and compare the experiment 
data with the results of the model developed in this paper will be a challenging work. 
Therefore in this paper the quark distribution is not discussed. However, the method 
developed in this paper can be easily extended to the evolution equation for the quark 
distribution. Even the method can also be applied to discuss the parton distribution at 
large-x, low region, because the DGLAP equation modified by the nonperturbative 
contribution can be applied to both the large-x region and small-x region (in small-x 
region, it is possible that the GLR equation describes the parton behavior , but it is 
difficult to check the GLR equation because the nonperturbative contribution weakens 



the shadowing effect showed in this paper). 



The main purpose of this paper is to developed a method, which the evolution equa- 
tion itself includes the nonperturbative contribution. The method has the important 
theoretical feature that it suggests a smooth transition from the low ("soft" ), where 
nonperturbative contribution dominates, to the large ("hard" ) region, where the per- 
turbative contribution dominates and the nonperturbative contribution can be neglected. 
Although the method is only a first step in considering the nonperturbative contribution 
to QCD dynamical equations, it can be checked not only by the HERA experiments but 
also by other strong interaction experiments because the method proposed in this paper 
is adaptable to wide region. 
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in part by the Foundation of Shanghai Jiaotong University. 

References 

[1] L.N.Lipatov, Phys. Rep. 286 (1997) 131 

[2] B.Badelek, J.Kwiecinski, Rev. Mod. Phys. 68 (1996) 445 

[3] B.Badelek, K.Charchula, M.Krawczyk, and J.Kwiecinski, Rev. Mod. Phys. 64(1992) 
927 

[4] L.V.Gribov, E.M.Levin and M.G.Ryskin, Phys. Rep. 100 (1983) 1 
[5] E.M.Levin and M.G. Ryskin, Phys. Rep. 189 (1990) 269 

[6] I.Abt et al, Nucl. Phys. B 407 (1993) 515; M.Derrick et al, Phys. Lett. B 316 (1993) 
412; 

[7] ZEUS Collaboration, M.Derrick et al., Phys. Lett. B 407 (1997) 432 



[8] M.Gliick, E.Reya and A.Vogt, Z. Phys. C 48 (1990) 471; C 53 (1992) 127; C 67 (1995) 
433 

[9] V.N.Gribov and L.N.Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438, 675; G.Altarclli and 
G.Parisi, Nucl. Phys. B 126 (1977) 298; Yu.L.Dokshitzer, Sov. Phys. JETP 46 (1977) 
641; L.N.Lipatov, Sov. J. Nucl Phys. 20 (1975) 93; 

[10] L.V. Gribov, E.M. Levin and M.G. Ryskin, NucL Phys. B 188 (1981) 555; A.H. 
Mueller and J. Qiu, Nucl. Phys. B268 (1986) 427 

[11] J.Kwiecinski, A.D. Martin and P.J. Sutton, Phys. Rev. D 44 (1991) 2640 

[12] J.Collins and J.Kwiecinski, Nucl. Phys. B 316 (1989) 307 

[13] K.J.Eskola, J.Qiu and X.N.Wang, Phys. Rev. Lett. 72 (1994) 36 

[14] D.V.Shirkov and LL.Solovtsov, Phys. Rev. Lett. 79 (1997) 1209 

[15] Yu.M. Shabelski and D. Treleani, Phys. Lett. B 403 (1997) 364 

[16] K. Golec-Biernat, M.W. Krasny and S. Riess, Phys.Lett.B 33 (1994) 

[17] A.L. Ayala, M.B. Gay Ducati and E.M. Levin, Nucl. Phys. B 511 (1998) 355 

Figure Captions 

Fig. 1. The dashed curve corresponds to the results of Eqs. (10) and the solid curve 
corresponds to the results of the DGLAP equation without considering the nonperturba- 
tive contribution. 



Fig. 2. As for Fig. 1, except for the GLR equation in the case R = 2eV^. 




Fig. 1, Bo He, Nonperturbative contribution to.. 
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